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ABSTRACT 


A veat-time processor fov NASA 16.3 single-polarized 
seatterometer data has been developed at the Remote Sensing Center f 
Texas ASM University. This processor is a hardware /software sys- 
tem which incorporates a microprocessor design and software for the 
calculation of normalized radar cross section in real time . Inter- 
face is provided to decommutate the NASA ADAS data stream for air- 
craft parameters used in processing and to provide output in the 
form of strip chart and pcm compatible data recording. The unit 
has been delivered to NASA/ JSC. 


AN AIRBORNE RADAR SCATTEROMETER 
SIGNAL PROCESSING SYSTEM 


INTRODUCTION 

The Remote Sensing Center at Texas A§M University 
has developed for the National Aeronautics and Space 
Administration (NASA) a real-time signal processor to be 
used in conjunction with the NASA/ Johnson Space Center 
(JSC) 13.3 GHz single polarized radar scatterometer system. 
It was the object of this development effort to provide 
a processing capability to speed the evaluation of scat- 
terometer data in earth resources investigations and provide 
data in a form that shortens the processing time required 
to produce a final product with respect to existing data 
reduction procedures. 

The effort was conducted in three phases. 

Phase I was a definition study and provided for the 
specification of the hardware design. Phase II provided 
for the fabrication of hardware and the development of 
test procedures and documentation. Phase III provided 
for the checkout and testing of the signal processor 
and for the delivery of the processor unit to NASA/JSC. 

This report constitutes the final report for 
the project. In this report the general requirements 
for the processing of scatterometer data are outlined 
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and the general theory of operation of the developed 
processor is developed. This report, together with the 
Processor Operation and Maintenance document, constitutes 
the bulk of the project documentation. 
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PROCESSOR SPECIFICATIONS 

The Real-Time Airborne Scatterometer Signal 
Processor (RASP) is a developmental extension of previous 
work done in real-time processing of scatterometer signals 
at Texas A§M University (TAMU) . This work was done for 
the Naval Surface Weapons Center (NSWC) (formerly Naval 
Ordnance Laboratory) and resulted in hardware designs 
and implementation of analog processing systems for the 
real-time classification of sea ice types from 13.3 GHz 
scatterometer measurements [1], This previous development 
and subsequent advances in large scale integration in 
solid state electronics, and microprocessor development 
established the feasibility of the design and construction 
of a real-time processor unit using combined analog and 
digital (hybrid) signal processing techniques. 

The project to design and construct the RASP 
unit was undertaken to bridge the gap between previous 
processor development and the availability of a fully 
developed processor system which could provide for in-flight 
evaluation of the techniques of real-time signal processing. 
This unit could a3.so provide a signal reduction capability 
for the 13.3 GHz scatterometer system so that processed 
results from data missions would be available for prompt 
evaluation of mission effectiveness in acquiring suitable 
scatterometer data. This evaluation of processed data 
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would also identify flight segments for ;*hich more precise 
data product could be obtained from other processing capa- 
bilities such as the existing NASA scatterometer data 
reduction software or the software package under develop- 
ment at Texas A§M. 

The RASP processor was designed and constructed 
as a single chassis unit using printed circuit wire wrap 
boards, large scale integrated circuits and hybrid elec- 
tronics, see Figure 1. The design centered around an 
Intel 8080 microprocessor system. The RASP unit uses 
the latest technology available on a production basis and 
incorporates system designs which are advanced state-of- 
the-art in scatterometer signal processing. 

The RASP processor was implemented to interface 
with several aircraft data systems. Aircraft flight 
parameters are input in either ADAS or NERDAS data formats 
from NASA flight systems and processor data output is 
available for strip chart recording. The output data are 
also provided by the processor in a Bi-phase Level (BI0-L) 
pem compatible format for recording. Documentation of 
these data reduction results provides a computer compatible 
format for further analysis and data tabulation. 

The principal data output product from the RASP 
processor is the tabulation of normalized radar cross 
section or scattering coefficient as a function of incident 
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Real-Time Airborne Scatterometer 
Signal Processor (RASP) 
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angle which is aligned for an illuminated target cell. 
These coefficients are calculated dynamically using 
instantaneous "values of the aircraft flight parameters 
and radar signals. This is an improvement on previous 
processing techniques which primarily utilize average 
aircraft flight parameters. The concept of real-time 
processing of scatterometer data also presents a novel 
approach to time-of-f light evaluation of mission success 
and timely consideration of processed data results. 
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SYSTEM OPERATION 

The RASP system is designed to interface with 
the NASA 13,3 GHz single polarized scatteromcter system. 

An understanding of this instrument and its signal outputs 
is essential to the design of the processor system. The 
NASA scatterometer is an active microwave sensor, a radar, 
which illuminates a target area with microwave energy 
and receives the return which is scattered from the 
target. Subsequently, the return signals are modified 
and recorded for signal reduction and processing. 

The Radar Equation 

The radar equation is used to describe the 
viewing of a target by a radar system. The following 
development of the radar equation follows that of Moore 
in [2, equations 9-66 through 9-69, p.419]. 

Consider the geometry of a radar remote sensor 
as illustrated in Figure 2. The power density at the 
target is that from an antenna of transmitter gain 
radiating a power P t upon a target at range R, thus, 

p G* 

r f 

Power density at target = j 

4ttR 

The poitfer received at the receiving antenna is the power 
density of the target backscatter across the effective 




Active Remote Sensor Configuration 


F i gure 2 . 
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aperture of the receiving antenna, that is, 

P t G t 1 

Power received - P - 5 - . a . ? A 

r 4'rrR* 4 itk r 

where a is the effective backscatter area or cross section 
of the target, and is the effective receiver aperture. 

If the radar system employs a duplexed antenna, 
i.e., if the receiver and transmitter utilize the same 
antenna, or if antennas of equal parameters are used for 
receiver and transmitter, then the relationship 



reduces the equation for received power to give one form 
of the radar equation: 



P t G 2 X 2 a 
(4 tt) 3 R 4 


The quantity in the radar equation of most 
interest to investigators is the radar cross section a, 
since it Is the only parameter relating information about 
the target. The cross section is generally expressed 
as an average over the cell area. A, and has been demon- 
strated [3, 4] to be a function of several parameters. 
They are, 


a° = cj/A = f(F, P, 0, <f>, S, e) 
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where 

c° = normalized target cross section 
F = frequency of incident energy 
P - polarization of incident energy 
0 = angle of incidence 
4* = target aspect angle 
S = surface roughness factor, and 
e = target complex dielectric constant 

Remote sensing applications rely on the cross 
section's dependence on the latter two factors, the surface 
roughness and dielectric constant of the target cell. 

These two factors are indicative of important surface 
parameters. For example, soil moisture content affects 
the dielectric constant [5], and wind speed over large 
bodies of water has been related to the surface roughness 
[4, 6, 7]. 

Radar Scatterometers 

An important class of radar systems is the con- 
tinuous wave (CW) scatterometer . A CW scatterometer 
irradiates the target with a continuous beam of energy 
and receives the return energy at a separate antenna. 

This return power may be compared to the transmitted 
power to determine the scattering cross section. The 
NASA 13.3 GHz scatterometer system, for which the data 
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processor was developed, is one of this class of scatterom 
eters. This system illuminates the area beneath the air- 
craft with an antenna beam plus and minus 60° from the 
nadir along track and 3° wide across track. Figure 3 
illustrates the area illuminated by the NASA system. 

As the name implies, CW Doppler systems utilize 
the Doppler effect; that is, the change in frequency of 
energy emitted from a source because of a velocity differ- 
ence between the source and the receiver. For an airborne 
CW Doppler system moving horizontally over the earth, 
the Doppler effect causes a return spectrum of frequencies 
from the illuminated area, with positive Doppler shifts 

r 

from the fore targets (approaching), and negative shifts 
from the aft targets (receding) . 

The amount of Doppler frequency shift, f^, 
introduced into the radar return from a particular point 
on the ground is related to the incident energy frequency, 
f i > by 



where 

v r » relative velocity of target toward the receiver, 
and 

g 

c = propagation velocity, 3 x 10 m/sec. 
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Since the relative velocity can be positive or negative, 
the return power is usually distributed in an asymmetric 
spectrum centered around the transmitted frequency. Thus, 
with a CW Doppler system, the return from a particular 
angle of incidence has a unique corresponding frequency 
in the return spectrum. 

Spectrum Folding and ’Sign-Sensing 1 

The processing of the data is difficult to 
achieve at the high transmitter frequency; therefore, 
the return signal is heterodyned at the radar receiver 
with a portion of the transmitter energy to produce a 
baseband signal. Unfortunately, this causes a "folding" 
of the fore and aft Doppler frequencies about zero Hertz 
in the frequency spectrum. In the NASA system, two base- 
band signals, termed the sine and cosine channels, are 
provided by splitting the received signal into quadrature 
signals before mixing. Figure 4 illustrates the nature 
of these two signals. Additionally, a signal proportional 
to the transmitter power is inserted into the cosine 
channel. Recovery of the fore and aft Doppler spectra 
(sign- sensing) and of the calibration signal is an 
essential step in the reduction of the 13.3 GHz scatterom- 
eter signals. 


Return Signal, S(.w) 



Sine Channel 




Figure 4. The NASA Scatterometer Output Signal 
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Spectrum Sampling and Cell Definition 

The recovered spectrum, although usually cen- 
tered around a lower and more convenient frequency, is 
otherwise identical to the original return spectrum. 

Each frequency corresponds to signal return from a unique 
angle of incidence. The return spectrum is sampled by 
a finite bandwidth filter which defines an area, or 
cell, on the ground. As shown by Jean [8], the sampled 
cell is outlined on the across track side by lines of 
constant Doppler shift having Doppler return frequencies 
equaling the upper and lower frequency limit of the band- 
pass sampled, while the along track sides are formed by 
the antenna beamwidth. The lines of constant frequency 
shift correspond to the intersection with the ground of 
cones of constant Doppler, concentric about the aircraft 
velocity vector. 

In sampling the return spectrum, an estimate 
of the power return from a cell sampled is obtained. 

This estimate is then used in the determination of the 
scattering cross section. By plotting these lines of 
constant Doppler shift, called isodops, for a constant 
frequency difference, Af , between each isodop, the cell 
definition and angular* dependence can be demonstrated. 
Such a plot is given in Figure 5. The isodops are seen 
to form hyperbolas which for a constant Af become more 
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widely separated as the angle of incidence increases. 

Thus, the cell length defined by the bandpass of the 
spectral sampling filter increases as the incident angle 
is increased. The cell area is then delineated by the 
intersection of the bandpass of the spectral filter and 
the limits imposed by the fan beam pattern of the antenna, 
A typical target cell is outlined in Figure 6. 

When the aircraft roll or drift is introduced, 
the area within the antenna pattern changes. In Figure 7 
is shown a typical target cell under aircraft roll or 
drift perturbations. The scattering cross section is 
calculated from this target cell. 


Coefficient Calculation 

The objective of scatterometer data reduction 
is to obtain the normalized scattering cross section, for 
a target cell, and ultimately to obtain a cross section 
function of a particular cell as a function of incident 
angle. Solution of the radar equation for c° identifies 
the parameters which are needed for the calculation of 
cross section. 


a°(e) 


A 2 G Z (0) ACO) P t 


The calculation of o°(0) involves known constants, system 
parameters, calculated geometric parameters and dynamically 



Figure 6. Instantaneous Ground Cell Area 
with no Aircraft Perturbation 
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Isodops 


Figure 7. Instantaneous Ground Cell Area 

with Aircraft Drift and/or Roll. 
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measured signal values. The system antenna gain parameter 
2 

G (9) is determined for the particular scatterometer system 
as a function of incident angle. The geometric parameters 
which appear in the solution of the radar equation are 
slant range to the target, R, and the illuminated cell 
area, A, These values are calculated using aircraft 
parameters which are dynamic rather than averaged values. 

The processor determines both the transmitted 
and received power by measuring the sign- sensed radar 
signals through a bandpass filter. For a given angle, 
a Doppler frequency based upon aircraft parameters is 
calculated and the bandpass filter is set for that fre- 
quency. The output of the filter is related to the 
received power P . The filters are also set for the 
12 KHz calibration voltage and its value is related to 
the transmitted power P^.. The ratio of powers required 
for the scattering coefficient calculation is obtained 



P r Ce)/-P t = 


71 


KE 


cal 


where 


E 

cal - measured voltage for a filter frequency of 
12 KHz 

E g - measured output voltage 



21 


H(£q) « the frequency response characteristic 
of the radar land/ sea filter, and; 

K = Ferrite Modulator Constant. 

Alignment 

The incident angle at which a target cell is 
viewed depends upon the relative position of the aircraft 
to that target. It is evident that the angle of Incidence 
to a fixed target changes with the flight of the aircraft. 
Thus, radar return at different viewing angles from a 
single target cell must be measured at different times. 
Consequently, aligning the calculated scattering coeffi- 
cients for a single target requires time shifting the 
calculated results. This time shift requirement is of 
major importance in providing a suitable data product 
from the scatterometer/processor system. 

Processing Requirements 

Thus, the processing requirements for the reduc- 
tion of 13.3 GHz radar scatterometer data are identified. 
The processor must: 

10 sign-sense the signals to unfold the data into 
fore and aft spectra; 

2) sample the power spectra at the appropriate 
frequencies to recover the radar return meas- 
urement and to obtain the measurement of the 
calibration level; 
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5) obtain the aircraft flight parameters so that 

geometric calculations of slant range and target 
area can be determined; 

4) calculate the normalized radar cross section 
from the solution of the radar equation; and 

5) time align the scattering coefficients from a 
single target so that the output from the pro- 
cessor gives aligned cross sections for each 
target area. 

This set of requirements defines, in general, the functions 
which must be implemented by the processor to accomplish 
the required data reduction and presentation. 

Functional Organization 

In order to take advantage of the speed and 
efficiency of analog sign-sensing and spectral sampling, 
along with the accuracy and flexibility of digital compu- 
tation, the RASP system combines both technologies to pro- 
vide a hybrid system capable of real-time calculation of 
the scattering cross section. In Figure 8 is shown the 
functional organization of the RASP. The analog processor 
accomplishes the sign- sense and spectral filtering opera- 
tions, while the digital processing section implements 
calculation, alignment and control functions. The digital 
processing section is developed around a microprocessor 
and, consequently, requires both hardware and software 
development. 
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HARDWARE 

The RASP system is designed to process signals 
generated by the NASA 13.3 GHz single polarized scatterom- 
eter. For real-time operation, it is necessary to accom- 
plish the required processing functions by implementation 
of a hybrid analog-digital processor system. Signals 
provided by the scatterometer system are preprocessed, 
measured and converted to digital form by the analog sub- 
systems. The digital subsystems implement a microprocessing 
environment which provides flexibility and precision in 
the computation of scattering cross sections. Additionally, 
the digital subsystem provides control for the selection 
of spectral sampling functions and digital conversion of 
signals at proper time intervals. The digital subsystems 
are developed to decode aircraft parameters available from 
a pem stream and to provide processing results in a pem 
stream compatible with NASA data reduction equipment. This 
section provides insight into the functional operation of 
the analog and digital subsystems. Additional detail is 
available in Processor Operation and Maintenance [3]. 

The Analog Subsystem 

The analog processor subsection performs the 
sign- sensing and spectral sampling functions. Spectrum 
recovery is achieved by remodulating in quadrature the 
twe baseband quadrature spectra output from the scatterometer 
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to a 20 KHz center frequency, then adding the two remodu- 
lated spectra. Figure 9 illustrates the total spectral 
separation procedure. 

A detailed block diagram of the "sign-sense ,r 
and AGC circuitry is shown in Figure 10. Inputs 1 and 2 
correspond to the sine (Redop 1) and cosine (Redop 2) 
channels of the scatterometer output. Unity gain buffer 
amplifiers provide the proper impedance termination for 
each channel. 

The two quadrature data channels are filtered 
to remove out-of“band noise by amplitude and phase matched 
active filters. Following the filters the signals are again 
buffered, and an adjustment is provided to remove any DC 
offset voltage introduced by the filter circuitry. 

The filtered sine and cosine channels are next 
mixed with separate 20 KHz signals which are themselves 
in phase quadrature. The resulting mixer outputs, both 
of which contain the folded fore and aft Doppler spectra 
centered at 20 KHz, are amplified and combined by sum and 
difference networks to yield the unfolded spectra. 

An AGC s; giial is derived from the sum and dif- 
ference channel outputs and is used to control the gain of 
the cosine channel such that the sine and cosine channels 
are maintained at equal amplitudes. 







sine 

channel 

in 


cosine 

channel 

in 
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The AGC voltage is produced by multiplying the 
sum and difference channels. This product yields a DC 
voltage proportional to the difference in the sine and 
cosine channel amplitudes plus higher frequency terms. 
This signal is filtered by a 20 Hz low-pass filter, 
amplified, and applied to the gain control input of the 
cosine channel post-mixer amplifier. 

The spectral sampling subassembly consists of 
a voltage controlled oscillator (VCO) , a balanced mixer, 
a 200 Hz low-pass filter, a linear integrator and an 
integrator as shown in Figure 11. In addition, there is 
a digital-to-analog converter which allows the digital 
processor to select the VCO frequency and an analog-to- 
digital converter which interfaces the output of the RMS 
module to the digital processor. 

The digital processor calculates the required 
Doppler frequency corresponding to the particular angle 
to be sampled within a given sample interval. The pro- 
cessor outputs a digital word to the DAC to provide the 
correct DC tuning voltage to the VCO. The VCO output 
is applied to the balanced mixer, shifting the desired 
portion of the Doppler spectrum to zero center frequency. 
The mixer output is provided with a DC offset adjustment 
to zero the DC output level for zero input signal. This 
output is then filtered by the 200 Hz low-pass filter. 



data 

channel 

input 


Balanced 



digital 

input 


Figure 11. 


Bloch diagram of spectral sampling and digital interface subassembly 
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The filter output, which is that portion of the processed 
signal from the desired Doppler resolution cell, is applied 
to the input of the linear detection and the output is 
integrated. 

After an integration period of 0.1 second, the 
signal is sampled by the ADC on command from the digital 
processor. The output of the ADC is a digital word which 
is proportional to the RMS voltage of the 200 Hz filter 
output. The square of this voltage is therefore pro- 
portional to the power received from the target within 
the Doppler resolution cell. 

The Digital Subsystem 

The primary function of the digital section of 
the RASP is the calculation of a 0 from the radar equation. 
This calculation, including computation of the cell area 
and range, is carried out in the Central Processing Unit 
(CPU). The CPU, actually a small digital computer, also 
provides control of the analog system in selecting the 
frequency setting of the VCO for sampling angle deter- 
mination and control of peripherals for data transfer such 
as a data recording device. 

Computer systems studied in the selection of the 
central processor ranged from minicomputers to micropro- 
cessors. Evaluation of the systems by Cariker [8] shows 


that minicomputers are actually more powerful than required 
and that the most desirable choice would be a microprocessor 
based CPU. Recent advances in integrated circuit technology 
have made possible these single chip processors with in- 
herently high reliability and economy. A study of the 
available microprocessors, with consideration of such 
factors as speed, computational power, availability, 
reliability, software and hardware support and simplicity 
of design, led to the choice of the Intel 8080 8-bit micro- 
processor as the base of the central processing unit. 

The Intel 8080 is an 8-bit central processing 
unit chip which has the capability of addressing up to 
64 K bytes of memory and up to 256 external devices. The 
8080 has a machine cycle execution time of 500 nanoseconds 
with three to five cycles required per instruction. 

The digital subsystem is developed on a bi- 
directional bus structure. Control lines are available 
to each device on the bus so that the CPU function can 
be affected by peripheral units. The CPU controls periph- 
eral units by way of an address bus and by certain control 
lines which indicate whether a read or write operation 
is being initiated. 

Various interface units to the CPU provide 
access to the data needed for the calculations. The 
analog subsystem functions as a peripheral t' - ’ the CPU 
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with, separate ports provided for the setup of the analog 
spectral filter center frequency and for the reading of 
the converted digital signals obtained from the analog- 
to-digital conversion unit. 

A separate digital controller and CPU peripheral 
accepts the aircraft state parameter from a serial digital 
data stream. The aircraft interface controller is capable 
of accepting either NERDAS or ADAS digital data. The air- 
craft interface continually updates its internal buffer 
memory with newly acquired data. During certain idle 
portions of the controller cycle, the buffer memory is 
available to the CPU to transfer the latest available 
aircraft parameters. Upon CPU initiated reads, the infor- 
mation is transferred to memory. 

There are several output modes for the processor 
system. Two serial digital interfaces are implemented to 
provide for data output. These include a standard com- 
munications interface which is compatible with Electronics 
Industry Associates (EIA) specification RS-232/C. This 
interface is secondary to the Bi-phase Level (Bi-0L) 
output used operationally. This particular format is 
compatible with NASA data acquisition systems so that 
processor output can be recorded under flight conditions 
and subsequently reduced and tabulated with standard re- 
duction equipment and procedures. Analog outputs axe 
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also provided for strip chart recording of processed 
scattering coefficients. 

Operator communication and programming is pro- 
vided by front panel functions. The front panel is 
divided into two sets of switches for operator inter- 
action, The top set of switches are primarily for system 
programming and diagnostic testing. These switches are 
normally unavailable to an operator during operational 
flight conditions. The lower half of the front panel is 
the operator's console. These switch functions and dis- 
plays provide the operator with a procedure for defining 
a particular processing configuration. The displays 
provide some monitor capability and warning lamps indicate 
failure of various parts of the subsystem. Operation 
modes are also controlled by an operation from the front 
panel. 
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SOFTWARE 

Processor software is implemented to provide 
control, computation and testing functions for the RASP* 
Various programmed modes are available and operator 
selection of these modes controls the processor functions. 
Computational sequences and data acquisition parameters 
are controlled by the processor software and preprocessed 
radar signals are converted to normalized scattering cross 
sections through software controlled calculations. Float- 
ing point arithmetic is accomplished through the soft\tfare 
package. Comprehensive diagnostic procedures are software 
controlled to assume proper operation of the RASP system. 

Modes 

The processor software operates in several basic 
modes according to operator selection. These modes of 
operation include Internal Calibration, Self -Test, Idle, 

Run, Halt and Reset. A flowchart involving these modes 
with brief descriptions is shown in Figure 12. The 
principal modes of interest are Idle, Self -Test and Run. 

The Reset mode allows the operator to reset the processor 
operation to a specific program initiation point. This 
reset is manually activated and can be accomplished during 
any phase of the processing. The Self-Test mode provides 
for a software guided procedure to test major hardware 
functions related to operator interface, memory malfunction. 








analog sign- sense, strip chart output and scatterometer 
calibration signal acquisition. The procedures for . 
operation of the Self-Test software are provided in 
Processor Operation and Maintenance [ 9 j - 

In the Idle mode, the operator defines the con- 
figuration setup for the processing to be accomplished. 

Angle selection, fore-aft target selection, and input and 
output channels are selected. This mode also allows for 
operator definition of aircraft flight parameters to be 
used for processing. Once a correct processing configura- 
tion is determined, the operator can initiate the Run 
mode for processing of data. 

When the Hun mode is selected, the processing 
of radar signals is initiated. The Run mode activates 
an acquisition cycle which is composed of a calibration 
step and ten coefficient calculation steps. Each step 
is accomplished in 100 ms; consequently, the acquisition 
cycle is 1.1 seconds in duration. During the calibration 
step, the filtered voltage corresponding to the scatterom- 
eter calibration signal at 12 KHz is measured. It is this 
value which is used for scattering coefficient calculation 
during the remaining cycle steps. It is also during this 
step that the bandpass filter is set up for viewing the 
signals at the next angle to be processed; the frequency 
of the filter for the next angle to be measured is calculated 
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the aircraft flight parameters for the calculation steps 
are determined; and the data obtained for the completed 
target cell are output. During the calculation phase of 
this step, signals are being filtered and integrated for 
the next angle to be processed. 

The coefficient calculation steps are identical 
except for the actual angle selected. The calculation 
step begins by sampling and analog-to-digital conversion 
of the filtered and integrated signals to be used in the 
scattering coefficient calculation. The integration 
time is fixed at 100 ms and the software waits until 
the integration is complete before the measurement is 
made. Once the signal is measured, the bandpass window 
of the processor is changed and integration at the next 
angle is initiated. While this integration is being 
accomplished, the data last acquired are processed. During 
this time, the frequency for the next angle to be filtered 
is determined; the measurement is converted to floating 
point; the power ratio is computed; target cell area is 
computed using aircraft flight parameters; the scattering 
coefficient is calculated; the coefficient is converted 
to decibels; and the data value is stored in an accumula- 
tion array for alignment. The step is thus completed 
and another calculation step is initiated until a total 
of ten calculation steps are taken. 
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In the tenth calculation step, the bandpass 
filter is set up to view the scatterometer calibration 
signal and upon the completion of this step, the calibrate 
step is initiated. 

The Run mode is active until an operator selects 
either Halt, Idle or Reset, at which time the Run mode 
terminates and the selected mode is initiated. The Halt 
mode provides for an orderly termination of processing 
steps where the Reset mode causes immediate termination. 

In the present implementation, there is little difference 
between terminating through the Halt mode or terminating 
directly into the Idle mode. 

A typical operator sequence involves turning 
the processor on and allowing a suitable warmup time. 

The next step is to select the Self-Test mode and to 
Reset the system. At this point, the Self -Test procedure 
is followed to assure that the processor is functional. 
Upon successful completion of the Self-Test procedure, 
the processor initiates the Idle mode. An operator 
selects the appropriate flight/processing configuration 
at this time and observes the processor confirmation 
of a correct setup. The next step is to initiate the 
Run mode at the beginning of the flight line. As the 
measurement line is completed, either Halt or Idle is 
selected to terminate the processing operation. With 



"9 


the processor in the Idle mode, the configuration can be 
modified to accommodate changes for the next flight line. 

Arithmetic 

The processor software is required to perform 
arithmetic and transcendental function over a broad range 
of numerical values. Calculation of the composite equation 
for the scattering coefficient requires numerous multiplies 
and divides. It also requires additions and subtractions 
but to a lesser degree. Floating point arithmetic is 
implemented for the processor due to the large dynamic 
range of the calculations. Since a larger number of 
multiplies and divides are required, the floating point 
representation is in a logarithmic form. Multiplication 
is accomplished simply by the addition of logarithms. 

Sixteen bits are used to represent floating 
point numbers. In the logarithmic form, the most signifi- 
cant bit of the number is the sign bit. When this bit 
is on, the number is negative. The next most significant 
seven bits are the characteristic of the logarithm in 
two's complement representation and the least significant 
eight bits represent the mantissa as an unsigned fraction 
with the binary point immediately preceding it. The base 
of the logarithm is two. 

A second form for floating point is used for 
the addition and subtraction of numbers. For the most 
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part, this form exists only internal to the addition 
and subtraction routines. This form is described here 
for completeness. The non- logarithmic form is similar 
to standard scientific notation where a number is repre- 
sented by a fraction and a power of the radix. In this 
representation, the 16-bit binary number has the most 
significant as a sign, the next most significant seven 
bits are the power of the Tadix two and the eight least 
significant bits are the binary fraction with the binary 
point assumed to precede them. The two representations 
are shown in Figure 13. 

Multiplication using logarithms is based on the 
relationship; 

log AB = log A + log B 
and division is accomplished by 

log A/B - log A - log B 

Multiplication is performed by binary addition of the 
logarithms of the floating point ; numbers which produces 
logarithms of the product. Likewise, division is deter- 
mined by finding the binary difference. Since the normal 
representation of the numbers is in a logarithmic form, 
multiplication and division software merely provides for 
the binary addition and subtraction. Results which would 
produce characteristics outside the maximum range of from 
+63 to -64 are set to the extreme value. The logarithmic 
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representation of zero is the smallest possible value, 
that is a logarithm with a characteristic of -64. Arith- 
metic functions and conversion tables return this value 
for zero when appropriate. 

It is noted that addition and subtraction are 
performed as true floating point numbers. Consequently, 
the number must be converted from the logarithmic form 
for the additions to be performed and returned to the 
logarithmic form for internal program consistency. The 
conversion between the two forms is accomplished by a 
table look-up procedure. The mantissa and binary fraction 
of the two forms are converted by the table and the 
characteristic and power-of- two are determined by the 
relationship, 

Power-of- two = characteristic + 1. 

Once the numbers are converted from logarithmic 
form, the fraction must be aligned according to the power- 
of-two. The smaller fraction must be logically shifted 
right and the power- of- two of the number incremented 
until the powers of both operands are equal. If the signs 
of the operands are different, the fraction of the negative 
operand must be complemented before addition. If the 
result of the operation is non-zero, and the most signifi- 
cant bit of the fraction of the result is not set, then. 
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post normalization is necessary. This is accomplished 
by shifting the fraction of an operand left and by decre- 
menting the power-of-two until a one appears in the most 
significant bit position of the fraction. At the end 
of the addition/subtraction routine the answer is converted 
to logarithmic form. 

Functions like sine, cosine and tangent are 
look-up functions where the input argument is expressed 
either in degrees or tenths of degrees in an 8-bit 
unsigned integer representation and the appropriate 
function routine is called. The entry corresponding 
to the value of the desired function is retrieved. A 
square root function is implemented which divides the 
logarithm of the floating point number by two. The 
routine returns a logarithmic zero for operands less 
than zero. 

The subroutines for the floating point arithmetic 
operations were initially written in Texas Instruments 
980A assembly language. Since the 980A supports FORTRAN, 
a test program was written to compare the results of the 
microprocessor simulation with known reliable results. 

Once the results from the two software packages agreed, 
the floating point software was converted to the micro- 
processor assembly language, retested, and integrated into 
the processing software. 


Area Calculation 


A critical parameter in the calculation of 
scattering cross section is the area of the target reso- 
lution cell viewed by the radar system. As previously 
discussed, the target cell is limited by the Doppler con- 
tours corresponding to the bandpass edges of the spectral 
sampling filter and by the limits imposed by the antenna 
pattern. These cell boundaries change as a function of 
aircraft altitude and velocity as well as aircraft atti- 
tude parameters of roll and drift. Pitch does not affect 
the area definition. 

The most exact method for calculating the target 
area is to perform a numerical integration over the cell 
boundaries; however, this is unreasonable in light of 
real-time processing constraints. When certain approxi- 
mations are applied to the problem, a closed-form solution 
for calculation of area is obtained. An approximate area 
is determined by obtaining the width of the terrain cell 
at the nominal look angle, as determined by the beamwidth 
of the antenna pattern, and multiplying this by the length 
of the cell as determined by the intersection of the Doppler 
contours associated with the band edges of the spectral 
sampling filter. Equations for each of these curves are 
developed in [8] and the resulting area term is defined as 
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In these equations, 

h = aircraft altitude 
v = aircraft velocity 

0 = look angle measured from aircraft nadir 

({> - aircraft drift 

ij> - aircraft roll 

fj = antenna across track beamwidth 
_ 1 

a = sin A (tan /tan <&) 

Af = bandwidth of spectrum sampling filter 
A - wavelength of transmitted energy. 
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Coefficient Alignment Procedure 

Alignment of the scattering coefficients from 
a single target is a major processor requirement. A 
processing algorithm has been developed which makes the 
alignment of computed coefficients a straightforward 
procedure. The cycle time required for the sampling of 
signals from all of the viewing angles defines the time 
between target cells which have been processed. The 
number of major cycles to be completed before a cell is 
overflown provides a convenient cell index. This index 
is used as the pointer to an array in which the values of 
the scattering coefficients and certain aircraft parameters 
are stored which relate to the indexed target cell. Cell 
alignment is accomplished by saving a computed scattering 
coefficient in the appropriate array position and reading 
out all the coefficients and aircraft parameters for a 
specific cell when all of its data have been acquired. 

In Figure 14 is shown a diagram of the cell 
acquisition sequence. Each ground cell is referenced by 
a sequence number which is the pointer into the array 
where values associated with that target cell are saved. 
This array is shown in Figure 15. Note that all of the 
scattering coefficient values for the fore angles 

o o 

(°0 " °3 ) have been acquired for the ith target cell. 




Figure 14, Cell Acquisition Sequence 
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Figure 15. Coefficient Alignment Array 
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During the next acquisition cycle, the aircraft 
advances along the ground track by one target cell, and 
the index i is increased by one. In other words, the 
nadir point of the aircraft is at cell i + 1. The radar 

o 

acquires cell i + 4 at the largest angle and for the 

i + 4th cell is computed. Similarly, for the i + 3rd 

o o 

target cell, 0 2 is computed, for the i + 2nd target, 

o 

is computed, and for the i + 1st cell is computed. 

These coefficient values are saved in the alignment array 

i 

corresponding with the index of the target cell and the 
particular angle associated with the coefficient. The 
array thus fills with each new acquisition and as the 
aircraft overflies the ground cell, associated flight 
parameters are also saved. 

At the time of o /erf light, all the information 
about the target cell has been acquired for the forward 
look direction. If this is the total information required, 
then the aligned scattering coefficients for this target 
cell can be output along with the associated flight param- 
eters at that point in time. Should the array position 
be continually reserved for a single cell, then as addi- 
tional cells are acquired, the array size grows without 
bound. However, after the data are output, there is no 
further need to save their position in the alignment 


array. 
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Modulo arithmetic is implemented in the alignment 
algorithm for the index pointer into the array. If the 
array is of a fixed size large enough to accommodate all 
target cells which are viewed by the radar at a single 
instant, and the pointer to the array is computed modulo 
that number or larger, then as a data set is output for 
a particular target cell, this array position is assigned 
to a newly acquired target cell. In this manner the array 
is of a fixed and, thus, bounded size, and an efficient 
utilization of storage memory is accomplished. 

The accommodation of aft data is a simple exten- 
sion of the procedure outlined and requires an array 
which is twice as large. The minimum size of the accumu- 
lation array for fore and aft data is given by the formula 

O 

XT 2h tan 60 
N = 

vT 


where 


N = the number of indexed vectors 
h = aircraft altitude in meters 
v = aircraft velocity in meters/second; and 
T = length of the acquisition cycle (sec) . 


For a fixed array size, the value of N is fixed and, con- 
sequently, each set of flight parameters must be checked 
to see if the array size is exceeded. The implemented 
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software checks this bound and if it is exceeded, defaults 
to a condition where the data are output as acquired 
without further alignment. 

Scattering coefficients which have an angle 
causing them to fall between the discrete target cell 
positions are assigned to the nearest target cell. Con- 
sequently, the average misalignment of scattering coeffi- 
cients is not greater than half the target cell separation 
distance. In Figure 16 Is shown a graph of permissible 
values of v and h for the specific software installation 
delivered with the procetjor. 
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CONCLUSIONS AND RECOMMENDATIONS 

A Real-Time Airborne Scatterometer Data Processor 
(RASP) has been developed to provide for reduction of the 
signals from the NASA 13.3 GHz, single polarized radar 
scatterometer system. This processor unit has extended 
previous techniques for the real-time hardware processing 
of scatterometer data. Previous development stopped short 
of actual scattering coefficient calculation. This develop- 
ment has provided for that calculation showing methods 
and procedures for its implementation. The unit has been 
developed beyond a feasibility demonstration model and 
actually provides for engineering evaluation of the 
processing techniques under limited flight conditions. 

Design goals for the system have been exceeded and the 
final unit provides improved flexibility, precision and 
data format. 

Other conclusions can be made. A design goal 
of 100 ms steps in signal measurement and the calculation 
of the scattering cross section has proven to be a con- 
servative number. It is estimated that an approximate 
15-20 pei cent improvement in the time performance of the 
processor algorithms could be met primarily by shortening 
the per point integration time. This improvement would 
result in the availability of closer target cells. 
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Observations of the total s/stem have identified 

o 

some problem areas. An unbalanced 60 Hz noise exists on 
one channel of the radar system. This signal is at times 
large enough to drive recorders into saturation for high 
noise peaks. Additionally, the existence of this noise 
on one channel and not on the second channel tends to 
defeat the AGC design implemented in the processor. 

Additionally, it has been observed that the VCO in the 
spectral sampling circuit is not as stable as would be 
desired. The resolution of the processor is reduced 
to make up for this uncertainty. 

Since this unit was involved in a development, 
it is more flexible than would be required for a strictly 
operational system. It is recommended that should similar 
processors be developed, their complexity be reduced from 
that of the current processor system. The programmer 
panel functions would more appropriately be placed into 
a test set rather than on the front panel. It seems 
reasonable that additional self-test steps could be added 
to facilitate processor alignment and repair. Specifically, 
steps to set upper and lower limits for the VCO would be 
useful . 

Determination of the radar calibration level 
through an alternate means offers some advantages concerned 
with VCO drift and some disadvantages regarding common 
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channel gains. It is recommended that alternate methods 
for the extraction of calibration signals be investigated 
for further developments of this variety. If unbalanced 
noise continues to exist in the 13.3 GHz scatterometer 
system, it is recommended that consideration be given 
to the disabling of the processor AGC provided channel 
gains can be maintained at like values. 

In future systems, additional consideration 
should be given to implementation with microprocessor 
circuits other than the Intel design. It is felt that 
the state-of-the-art in this area is changing at such 
a rate as to make current evaluation of products obsolete 
in a short time. Although the Intel system is most 
viable at this time, changes in technology may change 
this position. 
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